Pair correlations of scattered atoms from two colliding Bose-Einstein Condensates: 

Perturbative Approach. 
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I. INTRODUCTION 

The analog of correlated photon pair production pj] has 
recently been demonstrated using atoms. Both molecular 
dissociation and four wave mixing of deBroglie waves 
Q have shown correlation peaks. As in quantum optics, 
such atom pairs lend themselves to investigations into 
non-classical correlation phenomena such as entangle- 
ment of massive particles [4, 5, 6, 7]and spontaneous di- 
rectionality or superradiant effects 0, [§1 . From the point 
of view of the outgoing atoms, the underlying physics is 
very similar and thus theoretical descriptions should be 
applicable to both processes. The experiment using four 
wave mixing of metastable helium atoms in particular 
has yielded detailed information about the atomic pair 
correlations. Efforts to treat the experimental situations 
are therefore highly desirable. 

Theoretically, the description of condensate collisions 
in the spontaneous scattering regime requires a formula- 
tion that extends beyond the mean-field model [loj . 
In previous work on spherical Gaussian wave packets, 
within perturbative approach, we have given analytical 
formulas for the correlation functions [11, 12]. 

In this paper we extend our method to anisotropic con- 
densates to give an analytic description of the correla- 
tion properties of spontaneously emitted atom pairs in 
a geometry much closer to and in good agreement with 
the experiment [sl]. Numerical approaches using trun- 
cated Wigner method [H, [l^ ^iid positive-P method 
[Tsl . IT^ \u\ have also been used, in particular to give in- 
sight into the stimulation regime where bosonic enhance- 
ment comes to play. 

Here, we use the model of colliding condensates to ex- 
amine two types of correlations. First we shall focus on 
atom pairs originating from the same two body scatter- 
ing event. These consequently have nearly opposite mo- 
menta. Thus we analyze the opposite-momenta correla- 



tions of atom pairs. Second, we examine two body cor- 
relations between atoms scattered with nearly collinear 
momenta, a manifestation of the Hanbury Brown- Twiss 
(HBT) effect [lil, [H [H [ll] . In both cases, the demon- 
stration of a two particle correlation requires a measure- 
ment of the conditional probability of detecting a particle 
at position ri given that a particle was detected at r2. 
This probability is proportional to the second order cor- 
relation function G^^-'(ri,r2) of the field 6 of atoms, i.e. 



(^t(ri)^t(r2)5(r2)(5(ri)). 



We shall pay particular attention to correlations in mo- 
mentum space and compare these results with experi- 
mental data of A careful comparison of a numerical 
treatment based on the positive-P method [l3] with the 
experiment [H indicated reasonable agreement, but one 
of the limitations of the method, the short collision du- 
ration which could be simulated, left some unresolved 
questions. In particular, energy conservation is a less 
stringent constraint for short collision times, and thus 
one can wonder about the role this constraint plays in 
the experiment. The treatment given here is not sub- 
ject to this limitation and also agrees fairly well with 
the experiment for most of the experimentally accessi- 
ble observables. One observable quantity however, the 
averaged width of the collinear correlation function in a 
direction orthogonal to the symmetry axis, disagrees with 
the experiment and with Ref. T^. In our treatment, it 
is precisely the requirement of energy conservation that 
is at the origin of the difference. At the end of the paper 
we shall discuss possible explanations of this discrepancy. 

Let us first describe the experiment in which a collision 
of two Bose-Einstein condensates of metastable helium 
produces a cloud of scattered atoms. A condensate of ap- 
proximately 10^ He* atoms is created in a cigar-shaped 
magnetic trap with axial and radial trapping frequencies 
of LOz = Iny. 47 Hz and = 27rx 1150 Hz respectively. 
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FIG. 1: Velocity space representation of the pair production 
experiment. Raman pulses generate counter-propagating con- 
densates which collide and expand into disk-shaped clouds 
along the z-axis. Atoms scattered during the collision expand 
to form a spherical shell of correlated pairs. Note that the 
orientation of the axes in this article differs from Ref. 0] 

Three laser beams are used to transfer the atoms into 
two counter-propagating wave-packets by a Raman pro- 
cess, with a transfer efficiency of about 60%. As the 
wave-packets counter-propagate with a relative velocity 
of 2vrec = 18.4 cm/s, atoms from the two clouds collide 
via s-wave scattering, populating a spherical shell in mo- 
mentum space often referred to as the "halo" p^.[20l[2l|. 
In the experiment, about 5% of the atoms are scat- 
tered. In addition to splitting the condensate, the Ra- 
man transition transfers the atoms into an untrapped 
magnetic sub-state. The transferred atoms thus expand 
freely, falling onto a micro-channel plate (MCP) detector 
that allows the three-dimensional reconstruction of the 
position of single atoms with an estimated efficiency of 
10% [l^ll^. Knowing the positions of individual atoms, 
the initial momenta, and the second order momentum 
correlation function of the cloud of scattered particles 
can be computed. The precision of the measurement is 
limited by the finite resolution of the MCP. This factor 
will be taken into account in our comparison between the 
theoretical estimates and the experimental results. 



II. MODEL FOR SCATTERING 

To make the comparison, we introduce a simplified 
model for atom scattering during a collision of two Bose- 
Einstein condensate wave-packets. In this model we as- 
sume that two counter-propagating wave-packets consti- 
tute a classical undepleted source for the process of scat- 
tering. This concept is introduced in analogy to examples 
in quantum optics, where a strong coherent laser field is 
treated as a classical wave and its depletion is neglected 
[131 ■ We shall simplify the model further on. Since we 
assume that the two colliding condensates remain unde- 
pleted, the population of the 5 field of scattered atoms 



should be small, as compared to the number of atoms in 
the condensates. In such a regime, a Bogoliubov approx- 
imation is often used [1^ HI] , leading to linearized equa- 
tions of motion for the quantum fields. In our case, the S 
field of scattered atoms satisfies the Heisenberg equation 
(for details of the derivation, see [Til. IT^') 

thdtSir, t) = -^S{r, t) + 2#Q(r, OV'-Q(r, t)P{v, t). 

(1) 

Here '0±Q(r,t) is the c-number wave-function of the col- 
liding condensates with mean momentum per atom equal 
to ±TiQ. Moreover, the coupling constant g = ^'^^ ° is re- 
lated to the atomic mass m and s-wave scattering length 
a of He*. 

To permit analytic calculations, we model the conden- 
sate wave functions ip±Q{x,y, z,t) as Gaussians: 

xexp(-^(.^+,V^(.T^f), (2) 

where N is the total number of particles in both wave- 
packets. The radial (ar) and axial (cr^) width of the 
Gaussians are extracted from the initial condensate 
wave-function which is calculated numerically from 
the Gross-Pitaevski equation using an imaginary time 
method. In practice we fit J dv^ J dvy\'^ o{v)\'^ with a 
Gaussian function oc exp(— w^/x^) and then use = 
Ti/(m Xz)- We define ar similarly. Here, for simplicity, 
we neglect the spread of the condensates during the colli- 
sion. This assumption seems reasonable because most of 
the atom collisions take place before the two clouds have 
had time to expand. 

It is useful to change variables and rescale the field 
operator 

^t^t rM^r J^J(r,t)->5(r,t) 
"^cr^ at' 

which simplifies the equation of motion ([1]), i.e. 

(3) 

wherea=^,/3 = Qa„7-f7. 

The condensate density in momentum space then 
reads, 

|vI;o(k)|^ = ^fl^xp [~(3'ikl+jX)] (4) 

The three parameters a, (3 and 7 fully determine the dy- 
namics of the field of scattered atoms. For N — 10^ and 
hQ — TOWrec wc have a — 1053, (3 — 227, and 7 = 0.05. 

We also find Xz — 0.004 Vrcc, (Jz = 39 ^m, x,. — 
0.0870 Vice and <Tr = 2 /xm. The parameter a is a mea- 
sure of the strength of the interactions between particles. 



3 



As such, it governs the fraction of atoms scattered into 
the halo. As a consistency check, in Appendix A we give 
an ahernate estimate of a in the experiment using the 
observed fraction of scattered atoms. 

In Section [ml we derive an analytical expression for 
the second order correlation function in the perturbative 
regime. It is still an open question whether, for these pa- 
rameters, the perturbative approach applies. We tackle 
this issue after the evaluation of the G^^^ function is Sec- 
tion nil CI In Section IIVI we compare the perturbative 
results with the experimental data of ^SJ- 



III. DERIVATION OF G<^' IN PERTURBATIVE 
REGIME 

We shall begin the analytical calculations with a defi- 
nition of the Fourier transform of the b operator 

3/2 

<5(r,i) 



This particular form of Fourier transformation "incor- 
porates" the free evolution of the field. Substitution of 
Eq.® into Eq.© gives 



exp 



where A 



<5t(k',i). 



kr — k^Gx + kyGy, and is a unit 



vector in i direction. The above can be integrated for- 
mally, giving 



S(k,t)^Aj^ dre-'f^^e-^' J dk' e"^ 



exp 



2fl2 



■{K 



f3' 



Since in the Heisenberg picture the scattered field re- 
mains in its initial vacuum state and the evolution of 
the d field is linear, the second order correlation function 
G^^^(ki,k2) decomposes into 

G^^\kuk2;t) - (5t(ki,i),5t(k2,t)5(k2,t)<5(ki,i)) 

- G(i)(ki,ki;i) • G(i)(k2,k2;t) + |G(i)(ki, k2; 

+ |M(ki,k2;t)|' (6) 

where M(ki,k2;t) — {S{'ki,t)5{'k2,t)) is the anomalous 
density and G^^) (ki, k2; <) = ((5t(ki, t)(5(k2, t)) is the first 
order correlation function. Below we calculate these two 
functions in the lowest order and for a time t — oo be- 
cause all the measurements are made long after the colli- 
sion has finished. We expand (5 in a series of perturbative 
solutions. 



S{k,t = oo) = ^(5(')(k). 



where in the lowest order we get 

/•oo /> 

,5(i)(k) ^A dre-'P-e-^" / dk' e^('='+'=")^ (7) 
(-2:|!(k.-Hk;)2-^(fc, + fc;)2) J(")t(k'). 



exp 



A. Anomalous density: ki ~ k2 correlations 

The anomalous density in the first order is expressed 

by 

M(ki,k2) = (<5W(ki)<5(i)(k2)). 
Using Eq.([7|) we get 



M(ki,k2) =^exp( -^(ki,,+k2,.)2 



/32 

X exp ( — Y{ki,z + 



ki,zf^ dr exp (-i/3Ar - r^) 



where A = /3 (^1 - ^i±^^ . This gives 



M(ki,k2) = -i^^^exp ( -^(fci,, + fc2,.)' ) (8) 



exp 



This expression shows that the anomalous density de- 
scribes the correlations of atoms with opposite momenta. 
In other words, it is non-negligible only when ki ~ — k2. 
If this condition is not satisfied, the exponential func- 
tions drop quickly. Comparing this expression to Eq.(|4]), 
we find that the widths of the anomalous density have 
the same anisotropy and are two times larger than the 
condensate density. Moreover, this expression shows that 
this function is also non-negligible only for A < 1. As 
(3 is large, A ^ 1 only when fci ~ 1 and fc2 — 1- This 
requirement expresses the conservation of energy in the 
collision of two atoms. 



B. First order correlation function: ki ~ k2 
correlations 



In the lowest order we have 

G«(ki,k2)- (^«^(ki)<5(i)(k2)). 

Using Eq.® and ((5(o)(ki)(5(")t(k2)) = 5^^\ki - k2) we 
get 



4 



poo poo p 

G(i)(ki,k2) = \A\^ / dr / dr' exp [-r^ - t'^ + i/3(T - r')] / dk exp 
Jo ^0 J 



X exp 



4 



exp 



Under the assumptions that the foUowing three condi- 
tions are satisfied 
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« 1, 



1 



« 1, 



(9) 



where u= (ki-|-k2)/|ki-t-k2| and = (ki^r+k2,r)/|ki + 
k2| refers to the radial component of u. 

We show in appendix|B]that the atomic density is given 
by 



GW(k,k)^ „/i2:' , exp 



32V2n\ur 



2/3V(fc-l)^ 
|uJ2 



, (10) 



and the first order correlation function by 



G«(kr,k2) = ;;:i^exp 



32V27r|u 
X exp (^-^{uAkf^ ^1 - erf 

X exp ' 



Ak2 - ^Afc 



i/3uAk\ 

2V2 ; 



(11) 



We have introduced Ihi^ = 1 + aK, Ak ^ ki - k2 
and assumed |Ak| is small. 

The conditions ^ are fulfilled in the experiment of 
Ref 0] because the region ~ corresponds to the 
location of the two condensates and has been excluded 
from the analysis. The density of the scattered particles 
is peaked around k = 1 with a width of <C 1. We thus 
expect an anisotropic halo thickness, but the anisotropy 
is only strong around Mr- ~ 0, a direction which was in- 
accessible in the experiment of Ref. [1] 

As in the case of the anomalous density M, we can de- 
compose G^^-* (ki, k2) into factors expressing momentum 
conservation (1st line of Eq. 11) and energy conservation 
(2nd line of Eq. 11). We find that the widths of the mo- 
mentum contribution are V2 larger than the correspond- 
ing ones for M(ki, k2) [l3,[ll|- As discussed in Re fs. [15] 
and [16], the \/2 is due to the assumption of a Gaussian 
density profile. The energy contribution happens to be 
much more constraining than for M(ki,k2) because of 
the term uAk. If uAk = 0, meaning ki — k2, the width 
of G^^^(ki,k2) is given by the momentum contribution. 
But, if uAk 7^ and for instance if u is parallel to Ak, its 
width is oc 1/(3 even in the radial plane, in contradiction 
with the simple model developed in Ref. Q . 



C. Applicability of perturbation theory 



Perturbation theory is valid provided the scattering of 
atoms is spontaneous. When bosonic enhancement comes 
in to play, the perturbative approach fails. Here we give 
a simple estimate for parameters such as the number of 
scattered atoms and the dimensionless parameter f3 for 
which the perturbation is small and the above first order 
results can be used. 

A coherence volume can be attributed to each scattered 
atom. It is a volume in momentum space in which the 
atom is first-order coherent. In other words, if we choose 
a scattered atom with momentum k, the volume set by 
all the wave- vectors ki for which |G^^-'(k, ki)| is not neg- 
ligible is the coherence volume. If two bosons scatter in 
such a way that their coherence volumes overlap, their 
joint detection amplitude is enhanced by an interference 
effect. In other words, scattering into an already occu- 
pied mode is stimulated. The function G'-^-' permits an 
estimate of both the number of scattered atoms and their 
associated coherence volumes. If the number of scattered 
atoms is small, coherence volumes are unlikely to overlap, 
and stimulated scattering is negligible. In this situation 
we expect our perturbative solution to be valid. 

The above argument was used in the case of the colli- 
sion of two spherically symmetric (7 = 1) Gaussian wave- 
packets [121] and, in comparisons with numerical solutions 
of the equation for the field S proved to be correct. Here 
we apply an analogous reasoning for the case 7 7^ 1. A 
conservative estimate for the maximum number of scat- 
tered atoms for which the perturbative approach applies 
is A'crit = V/Vc, where ^ is a lower bound on the k- 
space volume into which atoms are scattered, and Vc is 
an upper bound on the coherence volume of an individual 
atom. 



In the comparison with the experiment (section IIV[) 
we analyze a fc-space volume which excludes angles 9 
smaller than tt/A. From Eq. pTI)) one sees that the density 
of scattered atoms is peaked around k = 1 with an rms 
width of sin 0/7/3. In the volume Q, the minimum rms 
width of the shell is (7/3\/2)^^. Taking twice this min- 
imum rms as the thickness of the shell, we find a lower 
limit on the volume of > Att/^(3. 

The analysis of Eq. (Ill|l shows that Vc reaches its 
maximum in Vl for ~ 7i"/4 (or ~ 37r/4, but due to 
symmetry we will focus on one of these values). If we set 
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7r/4 + SO, (f = 5ip and fci = ^2 = 1 we find 



This gives an angular area of coherence approximately 
equal to 87r/7/?^. Now we need to find the coherence 
width in the radial direction. Setting ki = (1 + (5fc/2)k/fc 
and k2 = (1 — (5fc/2)k/fc we get: 



exp 



The limit on the coherence volume is therefore 
647r/37/33. 

Combining the estimates of V and Vc, we find that 
critical number of atoms is given by A^crit = For /3 = 
227 we get A'crit ~ iC. In the experimental realization, 
the number of atoms detected in fi varied from 30 to 300. 
Assuming 10% detection efficiency this gives a maximum 
of 3000 scattered atoms. Thus the experiment should be 
in the perturbative regime. A similar argument is given 
in Ref . leading to a similar value of A^crit • 



IV. COMPARISON WITH EXPERIMENT 

The formulae ([5]) and pip cannot be directly com- 
pared with experimental data. This is due to an extra 
step which is made during the measurements: the joint 
probabilities measured in experiment are averaged over a 
region of interest Q which excludes the unscattered con- 
densates. We approximate ft by 9 £ [j, S [0, 27r] 
(where u = (sin 6 cos ip, sin 9 sin ip, cos 6) ) . 

In case of local momentum correlations, the normal- 
ization procedure is done by choosing ki and k2 almost 
equal: ki — k2 = (5k, where 5k is small. So we set 
ki = k -I- ^k/2 and k2 = k — (5k/2. The averaging corre- 
sponds to calculation of an integral 



(|G(i)(<5k)|2) = f dk |G(i)(ki,k2)p. (12) 

Then, this function is normalized by 

/ dk G«(ki,ki) •G(i)(k2,k2). (13) 

Let's denote the resulting normalized function by 
(|Sf^"'^-'(5k)p). As the anomalous density vanishes for local 
correlations, Eq.® gives 

5(^)(<5k) = l + (|.g«(5k)p). 

For (5k = we get g^^HO) = 2 [H]. 

In case of back-to-back momentum correlations, in 
analogy we have ki and k2 almost opposite: ki+k2 — Sk. 



We set ki = k -I- 5k/2 and k2 
the averaging corresponds to 



-k + 5k/2. Once again, 



(|Af(5k)|2) ^ f dk |M(ki,k2)|^ 



After normalization by function (|13p we obtain 
(|m((5k)p). For the opposite momentum correlations, 
G*-^-* vanishes, thus 

g(2)(^k) = l + (|m(5k)p). 

Let us now calculate the normalization function from 
(1131) . as it is common for both local- and opposite- mo- 
mentum correlations. From Eq. ljlip we have 



G(i)(ki^2) - 



327r\/2| Ul 2r 



exp 



1) 



2|ui, 



Now, in spherical coordinates, |ui,2r| = |sin0i^2|, where 
01,2 is an angle between vector ki,2 and axis z. Since ■^Sk 
is much smaller than k, we can approximate sin0i,2 — 
sin 6, where 6 is an angle between vector k and axis z 
and drop higher order terms in Sk in the exponentials. 
We end up with the approximate expression 



/ dk G«(ki,ki).G«(k2,k2)^-^ 



X I dk exp 
In 



- sm w 



li Sk ^ Sk ■ ex, k - Sk — kSk sin 61 cos <j) and if Sk — Sk ■ Bz, 
k ■ Sk = kSk cos 9. The resulting integrals are calculated 
numerically. 



A. Back to back momentum correlations 

As discussed above, we set ki = k + Sk/2 and k2 
-k + 5k/2. Using Eq.® 



|Af(ki,k2)P 



1 + crfi^ — 



2567r2 



exp 



2 2 ^ 4 



Here, A = /3(1 — k^ — ^). The averaging over fl is 
equivalent to 

(|A/(^k)H = 



dk |M(ki,k2) 



n 



2567r2 



xe 2 "'^^ 



dk e~T^ ( 1 + erfi2 ( ^ 



Numerical evaluation of this integral (for parameters j3 
and 7 as defined above) shows that the averaged anoma- 
lous density can be well-approximated by 



(|M((5k)|2) cx exp 



6 



As we see, the width of (|M((5k)p) is primarily deter- 
mined by the momentum conservation constraint, but 
the analysis shows that energy conservation plays a role, 
decreasing the predicted width in the xy-plane by of order 
10%. We normalize the second order correlation function 
by (|13|) and introduce an empirical parameter T^bb to ac- 
count for the fact that in the experimental data plots, 
the correlation functions are projections, and the fact 
that their heights were smaller than expected. We find 

This function is plotted in Fig. [21 using the value rjbh — 
0.032. We find good agreement with the experimental 
data in the x- and y- directions. In the z-direction, 
the width of the experimental peak is dominated by the 
detector resolution which is larger than the calculated 
width. 
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FIG. 2: (color online): Normalized opposite momentum cor- 
relations calculated in perturbative regime as compared with 
experimental data. Three plots correspond to three different 
directions. Here, Svi — h/m ■ Ski and Vrec = h/m ■ Q. 



radial integral. The radial one is 

ffc^dfcxexp^ 4/.V(^-l)^ 



The width of this Gaussian function is so small, that we 
can set k-dk ^ dk. Setting k — l + dk and extending the 
lower limit of the integral to -—oo gives Ir oc |ur|. Thus 



2ir 



3/4TV 



(|G(i)((5fc,)P) cxe-^*'^- / d^ I d0 

Jq Jtt/4 

( 0^ 

exp f — "^^^x ' "^(^i 



l + erfiM^-4tzi(e,¥') 



2V2 



where u{9,(p) = sin 6* cos (yS. This integral is calculated 
numerically and we obtain 

The result is again rescaled by the parameter rjd although 
it needs not to be identical to the back to back case: 

{g(^\Sk^)) = l + r^d{\g^^\6k,)\^). 

As = we deduce the value of rjci = 0.05. 

Now we set 6k — dkze^, and therefore u6k — 6kz cos 6. 
The radial integral is the same as in the previous case 
and we find 



(|G(i)(<5fc,)H ocexp 



3/4ir 



d0 



7r/4 



3'^ 

exp ( — —{Skz f ■ cos^ I 



1 + erfa I — cos 



2V2 



Numerically we find: 

{g^'HSkz)) = l + Vci{\9^'HSkz)\'). 



B. Local momentum correlations 

For the coUinear correlation function we choose ki = 
k+(5k/2 and ka = k~5k/2. Using Eq.((ni) and definition 
from Eq. (fT2ll we have 



(|G(i)((5k)|2 



111, 



dk 



n 2i^7r|u, 

■2 



■ exp 



(uSk)'' 



exp 



X exp 



Ski - ^Skl 
A(3^j^k^lf 



1 + erfi^ 



2V2 



uSk 



Let us now consider two separate cases. 

Let's set 6k — 5kxex. Then, u6k = 6k s'mO cos (p. In- 
tegration over the region consists of an angular and a 



We find that chosing ryd = 0.05 makes the observed 
heights match. 

Once again, because of the detector resolution, we find 
that the calculated peak is much narrower that the ob- 
served one in the z-direction. What is more surprising 
is that the widths of the correlation functions in the x- 
and y-directions are also narrower than those in the ex- 
periment. As can be seen from the discussion following 
Eg. pT]) . the peak width along the direction of the out- 
going atoms is strongly constrained by the energy con- 
servation requirement. This means that for scattering 
far from the z-axis (9 large) , the x- and y-components of 
the correlation function are narrower than they would be 
taking momentum conservation alone into account. This 
result contradicts the simple reasoning of Ref. 0] . In the 
next section we speculate about why neither the experi- 
ment nor the positive P simulation results reproduce the 
above width for the correlation function. 
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FIG. 3: (color online): Normalized coUinear correlations cal- 
culated in perturbative regime as compared with experimen- 
tal data. Three plots correspond to three different directions. 
Due to cylindrical symmetry of the colliding condensates, the- 
oretical results preserve this symmetry. Here, Svi = h/m ■ Ski 
and Vrec = h/m ■ Q. 



collisions take place before the clouds have had time to 
expand. Still, a quantitative estimate of the influence 
of the condensate expansion is another avenue for future 
analysis. 

Clarifying these questions may have ramifications be- 
yond atom optics. Conceptually similar experiments in- 
volving collisions between heavy ions have also uncov- 
ered discrepancies between observations and simple mod- 
els [13, [ill, the so-called "HBT puzzle". We hope that 
the work presented here will continue to stimulate care- 
ful thought about the four wave mixing process of matter 
waves. 
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V. CONCLUSIONS 

The perturbative result we have presented here, while 
rather complex, has the virtue that the results are ana- 
lytic and permit the identification the physical processes 
involved in the pair formation process. In particular the 
roles of energy and momentum conservation are clearly 
identified. Our results for the back to back correlation 
are in good agreement with the experiment. On the 
other hand the coUinear correlation function, as shown in 
Fig. [21 is in apparent contradiction with both the exper- 
iment and with the calculation of Ref . [l3| . The pertur- 
bative correlation function given in this work is narrower. 
This discrepancy clearly needs more attention, both the- 
oretical and experimental, but we wish to make some 
comments about possible causes. First, as discussed in 
Ref. flTl, the calculations using the positive P represen- 
tation are not able to simulate the entire duration of the 
collision; indeed only about 20% of the collision time can 
be simulated. Thus, energy conservation is not as strictly 
enforced leading to additional broadening in the calcula- 
tions of Ref. [l7|- Although this effect was discussed in 
that reference, the problem requires further scrutiny, it 
is not entirely clear to us which widths are most affected 
by a short collision time. Second, the experimental ob- 
servations are also subject to effects not treated here. It 
was briefly mentioned in Ref. Q that the mean field in- 
teraction between the escaping atoms and the remaining 
condensates may not be negligible. It is therefore im- 
portant to undertake an analysis of their effect on the 
correlation functions. Finally, an important simplifica- 
tion in the present treatment is the assumption that the 
condensates do not expand during the collision. This 
assumption seems reasonable because most of the atom 



APPENDIX A: DETERMINATION OF a 

When we introduced a, it was simply defined in terms 
of the number of atoms, the condensate size and the scat- 
tering length. Here we give a complementary estimate of 
a which provides a consistency check. The result es- 
sentially shows that our treatment is able to predict, to 
within experimental uncertainties, the number of scat- 
tered atoms. We start from Eq. pU)) . The integration of 
this equation over Q gives the number of scattered atoms 
to first order. This result, being a function of a, can 
be compared with the number of scattered atoms in the 
experiment. Knowing this number, we can evaluate a. 
First, using Eg. pTT)) . the number of scattered atoms in 
is given by 



M2 = 



32V27r|ur| 



dkexp 



2|uJ2 



Let us focus for a moment on the radial part of the above 
integral. 



^rad 



32\/27r u 



k dk exp 



(3^^\k^ - 1)2 



First, as the integrand is strongly peaked around fc = 1, 
the measured volume can be dropped, i.e. k^ ^ 1. Then, 
introducing k = 1 + 5k and assuming Sk is small we get 



-^rad 



32V27r|ur 



d{Sk) exp 



2|uJ2 



The lower limit can be extended to —ex:, giving 



'rad 



32V27r|u,. 



d{Sk) exp 



2\Ur 



a 7 
64 



8 



Integration over the angular variables gives a factor of of colliding Gaussians lies somewhere in between. This 



2^/2t: and 



TTa/S 2 2 

32 ' 



From the experimental data we know that the number of 
scattered atoms varies from 300 to 3000. For Nq, — 300 
we get a = 930 and for Nn = 3000 we get a = 2940. 
Thus the value of a = 1053 calculated from the model 



result, confirms that the choice of parameters such as 
and (Tr are reasonable. 



APPENDIX B: FIRST ORDER CORRELATION 
FUNCTION: ki ~ k2 CORRELATIONS 



To first order the G*^^^ function is. 



J 



nOQ nOQ n 

G(i)(ki,k2) = / dr / dr' exp [-T^ - t'^ + i/3(r - t')] / dk exp 

JQ Jo J 



{{ki^r+Kf + (k2,r +k^)^ 



X exp 



exp 



+ k^y ~ ^^{k' + kl)r 



Thus, in contrast to the anomalous density, we must lytically. Then, introducing x = and y = the 
perform a two-fold time as well as a three-dimensional grgt order correlation function is 
space integral. The space integral can be evaluated ana- 



GW(ki,k2) 



X exp 
X exp 

X 



167r3/2V2 

,,2 



exp 



|ki r — k2,r 



,2 



ki. 



k,J' 



-y + + + i 



4(1 + 2y2/(/?72)'; 
1 



kf + kl (ki+ks)^ 



exp 



da; / dy exp 

(ki,r + k2,r 



iV2y 



,2 y2 1 



exp 



iV2y- 



{ki,, + k2,zf y^ 



13 l + 2y2//32 







exp 





72(/372) l + 2y2/(/372)2_ 

_ (ki^^ + 'k2,rf ' 

4(l-F2y2//32) 



l + iV22//(/37^) ^i + ,^y/p 



where u = u,. + is a vector of unit length and di- 
rection ki + k2. As the scattering of atoms conserves 
energy and momentum, we expect that the density of 
atoms should be centered around |k| — 1 (which corre- 
sponds to |k| = Q in physical units). Moreover, from 



-y 



1 



M. 



the factor exp 

characteristic width of variable yisl/^/l-f 



we deduce that the 



Using the second of conditions 

,2 



we have 



exp 



-y 



1 + ^ 



u. 



exp 



u2' 

T . 



Since the characteristic range of y is 7/|ur|, all the terms 
proportional to y/ (3 and y/ 13^"^ can be dropped. This 
gives 



G«(kr,k2) = ,i^exp 



X exp 



167r3/2V2 



y_ 

7^ 



ki,z 

2 



^2,z| 



dx / dy exp 



2V2 



(kl,r 



^2\ 1 



(Bl) 



Now, by letting ki = k2 = k in Ea. (|Bl[) let us focus on 
the momentum density of scattered atoms. 



dy exp 



+ ipV2y {1 - k^) 
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X exp 



u^(l-fc^) 



X exp 



<(i-fc^) 



dy. 



From the above we deduce that the characteristic width 
of X is 1 which is much larger that the characteristic width 
of y. This allows another approximation - the limits of 
y integral can be expanded up from — oo to oo. The 
variables y and x effectively decouple, giving 



After integration over y and with k ^ 1, one obtains, 



G(i)(k,k) 



32V27r u 



■ exp 



2/3V(fc-l)' 



\\lr 



G(i)(k,k) = 



327rV2 



exp 



2 f 



Equation (jBl[) can be rewritten in the form 



G«(ki,k2) 



167r3/2y2 



exp 



1^ 
8 

■2 I e 



X exp 



L 7 



|ki r - k2,rP 



1^1, 



dx / dy exp 



2\/2 



u 



2 (kl,r+k2,r) 



Introducing ^'^^^'^^^ =1 + AK and Ak = ki — k2 , where 
lAkl is small we obtain 



g2 

X exp [ --^(uAk)2 



1 - erf 



i/3uAk' 
2V2 . 



G(i)(ki,k2) = 



a2/373 



32V27r|u 



■ exp 
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